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Abstract 

We examine three different ways a heavy quarkonium can dissociate at high 
temperatures. The heavy quarkonium can dissociate spontaneously when it 
becomes unbound at a temperature above its dissociation temperature. Fol- 
lowing the recent work of Digal, Petreczky, and Satz, we calculate the dis- 
sociation temperatures of heavy quarkonia taking into account the angular 
momentum selection rules and using a temperature-dependent potential in- 
ferred from lattice gauge calculations. We find that the selection rules change 
the dissociation temperatures substantially for charmonia but only slightly 
for bottomia. A quarkonium system in thermal equilibrium with the medium 
can dissociate by thermalization. The fraction of quarkonium lying above the 
dissociation threshold increases as temperature increases. A quarkonium can 
also dissociate by colliding with light hadrons. We evaluate the cross sections 
for the dissociation of J/^p and T in collision with vr as a function of the tem- 
perature of the hadron medium, using the quark-interchange model of Barnes 
and Swanson. We find that as the temperature increases, the threshold energy 
decreases and the dissociation cross section increases. 
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I. INTRODUCTION 



The suppression of heavy quarkonium production in a quark-gluon plasma has been a 
subject of intense interest since the pioneering work of Matsui and Satz Initial insight 
into the dissociation temperatures of heavy quarkonium was further provided by Karsch, 
Mehr, and Satz p[. Recently, Digal, Petreczky, and Satz reported theoretical results 
on the dissociation temperatures of heavy quarkonia in hadron and quark-gluon plasma 
phases. These are interesting results as they are related to the use of the suppression of 
heavy quarkonia production as a signal for the quark-gluon plasma 

The basic input of Digal et al. is the temperature dependence of the Q-Q interaction 
as inferred from lattice gauge calculations [^. We can understand conceptually such a 
temperature dependence of the Q-Q potential by placing Q and Q as an external source in 
a QCD medium. If the temperature of the medium is zero, the linear confining potential 
between Q and Q is the result of the alignment of the color electric fields at neighboring 
sites inside the flux tube, in analogy with the alignment of spins in a ferromagnet (Fig. la). 

(a) T = 



Q 




(b) High temperature 




(Fig. la) a schematic configuration of the color electric field be- 
tween Q and Q at T = 0. (Fig. 16) a possible configuration of the 
color electric field in a hot medium. Many different configurations 
of the color electric fields are possible, with a distribution centered 
about the mean color electric field configuration which is aligning 
along the flux tube. 

When the medium is at a high temperature, the gluon and light quark fields fluctuate 
and the alignment of the color electric fields due to the QCD interaction is reduced by the 
thermal motion for a random orientation of the color electric fields (Fig. 16). The greater 
the temperature, the stronger the tendency for thermal disorientation, and the weaker the 
resultant linear interaction between Q and Q. This tendency continues until deconfinement 
sets in at the phase transition temperature T^. 

The variation of the Q-Q potential gives rise to three different modes of heavy quarkonium 
dissociation at high temperatures: (1) spontaneous dissociation, (2) dissociation by thermal 
thermalization, and (3) dissociation by collision with particles in the medium. We shall 
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discuss them in turn. 

A heavy quarkonium dissociates spontaneously when the binding energy of the quarko- 
nium relative to a pair of final open charm or open bottom mesons vanishes. Utilizing the 
temperature-dependent potential for charmonia and bottomia, Digal et al. found that the 
dissociation temperatures for ip' , Xc, and T' in units of Tc are, respectively, 0.1 — 0.2, 0.74, 
and > 0.83, while the dissociation temperatures for J/ip, T, and Xb are above Tc [^3. 

Even at T = a proper description of the heavy quarkonium state should be based on a 
screening potential , as the heavy quarkonium becomes a pair of open charm or open 

bottom mesons when r becomes very large, due to the action of dynamical quark pairs. 
Previously, a screening potential was obtained to provide a good description of charmonium 
bound states, resonances, and e+e~ decay widths at T = pf. This screening potential can 
be used as the starting point for studying heavy quarkonia at high temperatures. Using a 
potential different from that of Digal et al. will help us find out the sensitivity of the values 
of dissociation temperatures on the potentials. We would also like to include spin- dependent 
interactions. Their inclusion is useful because they lead to a better description of the energy 
splittings of the quarkonium states which affect the dissociation temperatures. Furthermore, 
these spin and angular momentum quantum numbers give rise to selection rules for the final 
meson states and alter the dissociation threshold energies and the dissociation temperatures. 

The Q-Q potential at temperature T refers to the situation of placing an external quarko- 
nium (at rest) in a medium of temperature T. The changes of the gluon and light quark 
fields between the Q and Q due to the temperature of the medium give rise to a change of the 
Q-Q interaction and the subsequent change of the binding energy. The external quarkonium 
needs not be in thermal equilibrium with the medium. 

The quarkonium placed in a hot medium however will collide with particles in the medium 
and become thermalized. When the quarkonium reaches thermal equilibrium with the 
medium, the quarkonium will be in a mixed state whose occupation probabilities for the 
different energy levels will be distributed statistically according to the Bose-Einstein dis- 
tribution. There is a finite probability for this system to populate states whose energies 
lie above their respective dissociation thresholds. This fraction of the quarkonium system 
above the thresholds will dissociate into open charm or open bottom mesons. Dissociation 
of this type can be called dissociation by thermalization. Because the binding energies of 
the quarkonium states decrease as temperature increases, the probabihty of dissociation by 
thermalization also increases as the temperature increases. 

Finally, a heavy quarkonium can also be dissociated by collision with particles in the 
medium. The threshold energies and dissociation cross sections will be affected by the 
temperature dependence of the heavy quarkonium states. We would like to investigate how 
these dissociation cross sections depend on the temperature of the medium, using the quark- 
interchange model of Barnes and Swanson which has been applied successfully to many 
hadron-hadron [§-0 and hadron- (heavy quarkonia) reactions at T = [pr5|-p!8|. 

In Section II, we describe the Schrodinger equation and the potential used to calculate 
the energies and the wave functions of quarkonium states. We introduce a single-particle po- 
tential with simple spatial and temperature dependence to represent the results of the lattice 
gauge calculations. The single-particle states are calculated and exhibited in Section III. We 
discuss the selection rules for the spontaneous dissociation of heavy quarkonia in Section IV. 
In Section V, the dissociation temperatures are listed and compared with previous results 
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of Digal et al. [Q. In Section VI, we introduce the concept of dissociation by thermalization 
and estimate the fraction of the quarkonium system which can dissociate spontaneously as 
a function of temperature. In Section VII, the cross sections for the dissociation of ip and T 
in collision with pions are calculated as a function of temperature. We estimate the survival 
probability of a heavy quarkonium in collision with pions in Section VIII. We present our 
conclusions and discussions in Section IX. 

II. SCHRODINGER EQUATION FOR HEAVY QUARKONIUM STATES 

The energy e of the heavy quarkonium single-particle state {QQ)jls can be obtained by 
solving the eigenvalue of the Schrodinger equation 

|_V . -^V + V^i2(r, T) + A(r, T) j^j^^lr, T) = e{T)^jLs{r, T). (2.1) 

where 

A{r,T) = mi(r, T) + m2(r, T) - mi(oo,T) - 7712(00, T). (2.2) 

We have followed Ref. and describe the system as a two-body system whose properties 
and masses vary with r. At a distance r< R where R ~ 0.8 — 1.0 fm, the system consists of 
a quark and an antiquark, and their masses {mi(r, T), m2(r, T)} = {?7t,q (T), mQ(T)}. 

The nature of the the two-body system at r> R depends on the temperature. If the 
temperature is below the phase transition temperature Tc for light quark deconfinement, a 
heavy quark Q and antiquark Q separated at large distances will lead to the production of a 
light quark pair qq in between and the subsequent formation of open charm or open bottom 
mesons (Qq) and (qQ). At temperatures below the system at r> R therefore consists of 
a pair of heavy open flavor mesons, (Qq) and (qQ), and their masses {mi(r, T), m2(r, T)} = 
{MQg(T), MqQiT)}. By including dynamical quarks, the lattice gauge calculations of Karsch 
et al. incorporate this change of the configuration at large distances. In order to match 
with the lattice gauge calculations, the pair of mesons at large distances, (Qq) and (qQ), 
should be taken to be the lowest-mass pair, because the lattice gauge calculation results 
refer to those of the lowest-energy states of the system. The quantity fiu is the reduced 
mass. The interaction V12 is the interaction between Q and Q for R and between (Qq) 
and [qQ) for r> R. As the interaction between mesons has a short range, the interaction 
V12 (r,T) vanishes as r approaches infinity. 

Above Tc, light quarks are deconfined and the system is a separated Q and Q. Thus, 
{mi{r,T),m2{r,T)} = {mQ(T),mQ(T)} and A(r, T) = 0. The interaction Vi2(r, T) also 
vanishes as r approaches infinity because the interaction between Q and Q is screened at 
large distances in the deconfined phase. 

In Eq. ( p.l|) the energy e(T) is measured relative to the two-body pair at r — >• 00. 
The quarkonium is bound if e(T) is negative. The quarkonium is unbound and dissociates 
spontaneously into two particles if e(T) is positive, subject to selection rules which we shall 
discuss in Section IV. 

In the present manuscript, we shall limit our attention to T < for which the free 
energy contains contributions only from the color-singlet component and the extraction of 
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the color-singlet potential from the lattice gauge calculations is without ambiguity. For T 
slightly greater than T^, it is necessary to make assumptions on the relative fractions of the 
color-singlet and color-octet contributions in order to extract the Q-Q potential ||^ and the 
extracted potential there may depend on the assumed color-singlet and color-octet fractions. 

Limiting our attention io T < we note that the mass term A(r, T) can be approxi- 
mately represented as 

A(r,T) ^ [2mQ(T) - 2Mq,{T)] e{R - r). (2.3) 



The meson mass MqqiT) vary with temperature pO|-p^. Neglecting the kinetic energy of 



the heavy quark, the meson mass at temperature T is 

Mq,{T) = mgiT) + + [m,{T)f + (yQ,{T)) (2.4) 

which leads to 



[2mQiT) - 2Mq,(T)] - [2mQ(T = 0) - 2Mq,(T = 0)] = -2 [^p^ + [m,{T)f + {Vq,{T))\^^^ (2.5) 

In the evaluation of the Polyakov loop one obtains the free energy of the system from 
the lattice calculations. This free energy contains the interaction V12 between Q and Q 
or between {Qq) and (gQ), depending on the separation r. In addition, the free energy 
includes also the term on the right hand side of Eq. ( p.5| ) due to the variation of light 
quark q with temperature, since the light quark fields are dynamical variables and they 
changes the energy of the light quark, and the interaction between the light quark and the 
heavy quark {VqqiT)). The free energy of the lattice gauge calculations actually yields an 
effective interaction V{r,T) which includes Vi2{r,T) and the temperature dependence of 
A(r,T) - A(r,0): 

V{r, T) = \/i2(r, T) + [A(r, T) - A(r, T = 0)]. (2.6) 

This effective Q-Q interaction V{r, T) vanishes at large distances. 

In terms of this potential V{r, T) extracted from the lattice gauge calculation, the 
Schrodinger Eq. ( |2.1D is then 

|_V . -^V + V{r, T) + A(r, T = Q)\ijjLs{r, T) = e{T)ijjLs{r, T). (2.7) 

I 2/ii2 J 
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Fig. 2. Quark-antiquark potential as a function of r obtained by 
using Eqs. (pD-fCT) and bo = 0.35 GeV^. 

Our first task is to obtain a convenient representation of V{r, T) from the results of lattice 
gauge calculations of Karscli et al. 0. Following Karsch et al. 0, we represent V{r,T) by 
a Yukawa plus an exponential potential 0| 



y(r,T) 



3 r 



(2.8) 



where b{T) is the effective string-tension coefficient, /i(T) is the effective screening parameter, 
and the potential is calibrated to vanish as r approaches infinity. The results of the lattice 
calculations of Ref. |^ for T < Tc, which has been normalized to the Cornell potential 
at short distances, can be described by 



6(T)=6o[1-(T/T,)XT,-T), 



(2.9) 



and 



fi{T) = fioO{T,-T), 



(2.10) 



where bo = 0.35 GeV^, /io = 0.28 GeV, and 9 is the step function. The value of fio has been 
fixed to be the same as in T = 0, and the value of bo is close to the value of 6 = 0.335 GeV^ 
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obtained earlier at T = 0]. 
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Fig. 3. The potential difference between r^/a = 5 and r^/o^ = 0.5 
as a function of the temperature. The curve gives the results of 
Eqs. ( p.8|) and (|2.10|) with the parameter Bq = 0.35 GeV^ and the 
solid points are lattice gauge results of [Q. 

As a comparison, we plot in Figure 2 the quantity [V{ry/a) — V{0.5)]/ ^/o^ calculated 
with the above potential, Eqs. ( p.8|) - (|2.10| ), as a function of ry/a for charmonium at different 
temperatures. The scale constant ^/a■ = 0.425 GeV is used to convert distance and energy to 
dimensionless numbers as in the lattice gauge calculations of Ref. p. The general features 
of the spatial and temperature dependence of the potential in Eqs. (|2.8| )- (|2.10| ) agree well 
with the lattice gauge results [Fig. 6 of 0]. In particular, there is good agreement between 
the variation of the potential difference between r^/a = 5 and Vy/a = 0.5 as a function of 
the temperature calculated for charmonium, as shown in Fig. 3. 

In our calculation, we use a running coupling constant obtained in another meson spec- 
trum study [|TB|, 



127r 



{33 - 2n f) \n{ A + QyB^)' 



(2.11) 



where A = 10 and B = 0.31 GeV, with Q identified as the mass of the meson. The above 
formula gives a.^ ~ 0.32 for charmonium and as ~ 0.24 for bottomium. We include, in 
addition, the spin-spin, spin-orbit, and tensor interactions. 



spm— spin 



spin— orbit 



4 

- X — 

3 3mim2 



vr 



3/2 



+ ClsL-S + CtSi2, (2.12) 
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where Sj are the spins of the interacting constituents, L is the orbital angular momentum, 
S = S1 + S2, and S12 = {3{(Ti-r){(T2-r) — cri-cr2}. For our calculations, we take {Cls, Ct} = 
{34.6,9.78} MeV for charmonium, {Cls,Ct} = {14.25,2.98} MeV for bottomium. The 
width parameter d is the range of the spin-spin interaction. Following Godfrey and Isgur 
23| , we use a running width parameter to include the relativistic effects on the width of the 



spin-spin interaction, 

d' = ^o{i; + ^{ ^S'rZ\. ] K^nf^l. (2-13) 




where we obtain the values of ctq = 0.591 GeV and (Xi = 1.967 by examining the hyperfine 
splittings of n-p and rjc-J/ip. 

For T = 0, Mqq are known from experimental B and D meson masses. Using the 
V{r, T = 0) extrapolated from lattice gauge calculations, we vary rric and mj, and find the 
zero-temperature constituent quark masses = 1.89 GeV and rrih = 5.22 GeV which give a 
good fit to the heavy quarkonium spectrum (Figs. 4 and 5). These heavy quark constituent 
mass values are close to those of rric = 1.84 GeV and rub = 5.18 GeV in the Cornell potential 
I9| , whose short-distance values have been used to normalize the results of the lattice gauge 



calculations g]. 

The string tension term in the screening potential Eq. ( |2.8| ) needs to have the prop- 
erty that it vanishes at r — > 00, and the exponential potential has been conveniently cho- 
sen to describe such a property. It is of interest to compare the exponential potential 
—bo exp{— /ior}//io with the standard linear potential or at T = 0. The effective local string 
tension, as obtained by taking dV/dr, is a for the linear potential, independent of r. For the 
exponential potential, it is bQexp{—fior} which is equal to 60 at r = and zero at r — > 00. 
At the radius of J/ip, the local string tension is 0.57 60 = 0.20 GeV^, which is close to the 
standard value of a = 0.18 GeV^, as it should be. It is therefore reasonable that the pa- 



rameter bo is approximately two times the standard value of a |T8[. The large value of 60 in 
comparison with a was noted earlier by Ding et al. [0. 



III. HEAVY QUARKONIUM SINGLE-PARTICLE STATES 

The eigenvalues of the Hamiltonian can be obtained by matrix diagonalization using a 
set of nonorthogonal Gaussian basis states with different widths, as described in Ref. |p!8[ . 



For the construction of the Hamiltonian matrix, the matrix elements of the Yukawa potential 
and exponential potential between Gaussian basis states are 

^ /2^y- e.p{-^}. (2, . 5, , (3.1, 

and 
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where U is the cyhnder parabohc function, (5 is the width parameter of the wave function 
basis, and we have used the notation and normahzation of Ref. [18|. The matrix elements 
of the kinetic energy operator are 



where 



2/ii2(r; 



■v|j) = T,,> + r,,-< 



I 



-{21 + 1)2' Ii{z) + 2'+' h+,{z) 



(2Z + 1)!! 

tj (2/ + 3)!!-2'+2/,+2(^) 



(2/ + 1)!! 



(3.3) 



(3.4) 



T 



«+3/2 




(2/+ 1)!! 



(2/ + l)2'/K^)-2'+^/,+i(z) 



2'+2/z+2(^) 



(z+jT (2/ + 1)!! 



(3.5) 



li + j 



(3R, 



(3.6) 



IJz) 



IT JO 




1/2, 



A=l 



(3.7) 



and //> and /i< are the reduced masses in the region r > R and r < R respectively. 

After the Hamiltonian matrix is constructed using these matrix elements, the single- 
particle state energies and wave functions can be obtained by matrix diagonalization. 

In Figs. 4 and 5, we show charmonium and bottomium single-particle states as a function 
of temperature. On the left panel of each figure, the experimental single-particle energies 
at T=0 are also shown for comparison. The single-particle energies {ej} rise with increasing 
temperatures. At a certain temperature, the energy of a single-particle state will rise above 
its threshold for dissociation. But what is its dissociation threshold? 



IV. SELECTION RULES FOR SPONTANEOUS DISSOCIATION OF HEAVY 

QUARKONIUM 

In considering the dissociation below Tc, it is necessary to find the selection rules for 
the spontaneous dissociation of a heavy quarkonium state with initial quantum numbers J, 
Li, and Si into two mesons with a total spin S and a relative orbital angular momentum 
L. We shall only discuss the selection rules for charmonium states, as similar rules apply in 
the case of bottomium states. We limit our attention to final states consisting of a D or D* 
with an antiparticle D or D*. The parity of a quarkonium state is (—1)'^*+^. On the other 
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hand, the parity of a pair of dissociated open charm mesons is (—1)^. Parity conservation 
requires AL = |L — Lj| = 1 in this dissociation. Hence, J /if) and if)' will dissociate into a 
pair of open charm mesons with L = 1, while x states will dissociate into a pair of open 
charm mesons with L = or 2. 




Fig. 4. Charmonium single-particle states as a function of temper- 
ature. The threshold energies are indicated as horizontal lines. The 
solid circles indicate the locations of the dissociation temperatures 
for states with significant decay branching fractions into J/ip. 

The threshold for spontaneous dissociation of J/ip and ip' is therefore M{D) + M{D) + 
2/2/i£)/)i?^ where fi^i) is the reduced mass of the final two- meson system and the term 
2/2fij^[)R'^ represents the L = 1 centrifugal barrier height at the distance R at which D 
and D become on the mass-shell. The reduced mass fiDfj is equal to mo/^,. For numerical 
purposes, we use R = 0.8 fm. 

Consider now the dissociation of Xc2- If the final orbital angular momentum of the two- 
meson state is L = 0, then because J = L + S = 2, the final two-meson system must 
have S = 2 and the two mesons must each have S = 1. The threshold energy for the 
dissociation of Xc2 into two mesons in the L = state is M{D*) + M{D*). On the other 
hand, if the final state has L = 2, then any spin combination of the final meson is allowed, 
and the lowest threshold for the dissociation of Xc2 into two mesons in the L = 2 state is 
M{D) + M{D) + Q/2i^ddR^- 
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Consider next the dissociation of xa- If the final state has L = 0, then because J = 
L + S = 1, the final meson states must have S = 1. The threshold energy for the dissociation 
of Xci into two mesons in the L = state is M{D) + M{D*). The dissociation of Xi into 
D + D is not allowed. 




Fig. 5. Bottomium single-particle energies as a function of tem- 
perature. The solid circles indicate the locations of the dissociation 
temperatures for states with significant decay branching fractions 
into T. 

One can obtain the selection rules for rjc, Xco, and he in a similar way. We list the 
selection rules for the S'-wave charmonium states {J/ip, and rjc) in Table I, and the 
P-wave charmonium states {xcj and he) in Table II. 



Table I. Selection rules for the dissociation of S-wave charmonia. 



Initial Heavy Charmonium 


Final Mesons 


L 


Threshold 




D + D 


1 


M{D) + M{D) + 2/2fiDDR^ 




D + D* 


1 


M{D) + M{D*) + 2/2^iDD*R' 




D* + D* 


1 


M{D*) + M{D*) + 2/2/i^*5*^' 


Vc 


D + D 


1 


not allowed 




D + D* 


1 


M{D) + M{D*) + 2/2/i^5.i?^ 




D* + D* 


1 


M{D*) + M{D*) + 2/2fiD'5,R' 
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Table II. Selection rules for the dissociation of P-wave charmonia. 



Initial Heavy Cnarmonium 


final Mesons 


L 


i nresnold 


Xc2 


D + D 





not allowed 






2 


M[D) + M(L)) + 6/2ij,d£)R'' 




D + D* 





not allowed 






2 


M{D) + M{D*) + 6/2/i£,5*i?^ 




D + D 


U 








2 


M[1J ) + M[D ) + b/ IjJi H 


Xcl 


D + D 





not allowed 






2 


not allowed 




D + D 





M[U) + M[U ) 






2 


M(i:)) + M{D*) + Q/2hjjq*R^ 




D + D 





M[D ) + M(iJ ) 






2 


M[D ) + M[D ) + 6/2//£)*5*it^ 


XcO 


D + D 





M[D) + M[D) 






2 


X n J 

not allowed 




+ 





j_ 1 1 1 

not allowed 






2 


X 11 J 

not allowed 




L) + L) 





M[D ) + M[D ) 






2 


^(L' j + ^(17 ) + b/ 2iJ,j)*£,*K 


he 


D + D 





not allowed 






2 


not allowed 




D + D* 





M(D) + Af (L)*) 






2 


.U(D) + .U(D*) + (3/2/,,,.,,./?^ 




D* + D* 





M(L)*) + M{D*) 






2 


M{D*) + M{D*) + 6/2/XB.5*i?^ 



Similar selection rules can be written down for bottomium dissociation by replacing J/tJj 
by T, ijj' by T', Xc by Xb, he by hb, and D hy B in Tables I and II. 

V. DISSOCIATION TEMPERATURES FOR CHARMONIA AND BOTTOMIA 

Fig. 4 shows the charmonium single-particle states as a function of temperature. The 
dissociation temperatures of different quarkonium states can be determined by plotting 
the threshold quantities and the state energies of the quarkonia. The points of intercept, 
as indicated by sohd circles in Fig. 4, give the positions of the dissociation temperatures. 
Because the observed J/ip receives little feeding from the Xco and he states, and rjc is not 
observed in dilepton decays, we shall not include the dissociation temperatures of Xco, he, 
and r]c (and similarly Xbo, hb, and r]b) in our consideration. Fig. 5 gives similar temperature 
dependence of the bottomium single-particle states as a function of temperature. 

Table III. The dissociation temperatures in units of Te for various heavy quarkonia 
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TT 

Heavy 
Quarkonium 




Xc2 


Xcl 




1 


/ 

Af)2 




1 


Xb2 


Xbi 


1 


rji Irji 

J-dl J-c 


U.oU 


o.yi 


o.yo 


o.yy 


0.0 1 


n oo 


n oo 


n n^^ 

o.y6 


\ 1 nn 
> 1.00 


\ inn 
> 1.00 


\ 1 nn 




0.37 


0.66 


0.72 


0.94 


0.54 


0.75 


0.78 


0.95 


1.00 


>1.00 


>1.00 


(Selection Rules 
not invoked) 


























0.1-0.2 


0.74 


0.74 


1.10 


0.75 


0.83 


0.83 


1.10 


1.13 


1.13 


2.31 


(Digal et al.) 

























We list the dissociation temperatures of charmonia and bottomia in Table 111. If the 
dissociation threshold is just M{D)+M{D), the dissociation temperatures in units of for 
i^'j Xc2, Xcl, and J/ip in units of Tc would be 0.37, 0.66, 0.72, and 0.94, as hsted in the second 
row of Table 111. This is in rough agreement with the results of Digal et al. which give the 
values of 0.1-0.2, 0.74, and 1.1 for ip', Xc, and J/ip respectively. One of the main reasons 
for the observed differences is that the locations of the dissociation temperatures depend 
on the energies of the theoretical single-particle states. Hyperfine interaction, spin-orbit 
interaction and the shape of the potential affect the positions of the single-particle states. 
Although the fine structures are small in terms of the separation between major shells, they 
are nonetheless important in determining the locations of the dissociation temperatures. 

With the additional threshold energies due to the angular momentum selection rules, 
the dissociation temperatures in units of for Xc2, Xci, and J/ip are shifted to 0.50, 
0.91, 0.90, and 0.99 , as listed in the first row of Table 111. The increase in the dissociation 
temperature is greatest for Xc2 and least for J/ip. 

For the bottomium system, because the mass of the bottom quark is very large, the 
mass difference between B* and B and the centrifugal barrier energies are small; the selec- 
tion rules do not lead to substantial changes of the thresholds. The shifts of the dissociation 
temperatures are not as large as in the case of charmonium states. The dissociation temper- 
atures for various bottomium states are listed in Table 111. We confirm the general features 
of the results of Digal et al. but there are also differences in the details as the dissociation 
temperatures depend on the potential and interactions, as well as on the selection rules. 



With the temperature-dependent potential of Eqs. (|2.8|) - (|2.12| ), all the dissociation tem- 
peratures of heavy quarkonia are below T^, except for Xbi, Xb2, and T. The present analysis 
uses a potential valid in the range T < and cannot be used to determine dissociation 
temperatures beyond Tc. 



VI. DISSOCIATION OF QUARKONIUM BY THERMALIZATION 

Besides spontaneous dissociation, dissociation by interaction with particles in the 
medium can also lead to dissociation. In order to bring out the salient features, it is useful 
to present simplifying descriptions of these dissociation processes in terms of dissociation 
by thermalization and dissociation by collisions. A full description will involve the interplay 
between the dynamics of thermalization and the dynamics of different types of dissociation. 

In the dissociation by thermalization, we consider a two-step process. In the first step, 
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a thermal equilibrium can arise from inelastic reactions of the type 



h+{QQ)jLs^h' + {QQ)j>L>s' (6.1) 

in which the collision of a hadron h with a quarkonium in state JLS results in the exci- 
tation or de-excitation of the quarkonium state. If the heavy quarkonium reaches thermal 
equilibrium with the medium, the occupation probabilities of the heavy quarkonium state 
Ci will be distributed according to the Bose-Einstein distribution, 

= exp{(e.-/i)/T}-l- ^^-^^ 
If there is no dissociation of the quarkonium, then the sum of the occupation numbers is 

n = J2i2Ji + l)ni = 1, (6.3) 

i 

where 2Ji + 1 is the spin degeneracy. This equation can be used to determine the chemical 
potential fi. At temperature T, the fraction of heavy quarkonium whose single-particle state 
energies lie above their dissociation thresholds is 



/ = E 



(6.4) 



In the second step, after the quarkonium reaches thermal equilibrium with the medium, 
there is a finite probability for the quarkonium system to be found in excited states. If these 
excited states lie above their thresholds for spontaneous dissociation, the heavy quarkonium 
system will have a finite probability to dissociate into open charm or open bottom mesons. 

Previously, dissociation by thermalization was studied by Kharzeev, McLerran, and Satz 
2^ with a temperature-independent dissociation threshold and free-gas continuum states 



in a Boltzmann distribution. They found that the dissociation by thermalization does not 
provide a significant amount of J/ip dissociation at temperature T=0.2 GeV. Results in 
Sections JTT| and ^ indicate however that the dissociation thresholds and the positions of 
the single-particle states change with temperature. Furthermore, the set of states included 
in Eqs. ( |6.2| ) and ( |6.3| ) in the quarkonium thermal equilibration process ( |6.1[ ) should consist 
only of bound and resonance states [^, because the thermal equilibration of the heavy 
quarkonium QQ takes place when Q and Q remain in the vicinity of each other. Free-gas 
continuum states off the resonance represent Qq and qQ states whose spatial amplitudes lie 
predominantly outside the quarkonium system and they do not participate significantly in 



the quarkonium thermal equilibration process (6.1). 

To study dissociation of heavy quarkonium by thermalization, we evaluate the fraction 
/ of heavy quarkonium whose single-particle energies lie above their dissociation thresholds, 
as a function of temperature, using the single-particle states obtained in Sections |ITl| and 
[rV| . The results are shown in Fig. 6. The fraction / increases with temperature, due 
predominantly to the change of the positions of the single-particle energies. In Fig. 6, a 
state label shown along the curves denotes the onset of spontaneous dissociation of that 
state. As the temperature approaches T^, the fraction lying above dissociation thresholds is 
quite large for charmonium, and is considerably smaller for bottomium. 
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By the dissociation of the excited heavy quarkonium, the occupation number is no longer 
constant and the rate of dissociation of the quarkonium is 



^ = -^A.(2J, + lK 



(6.5) 



where Aj is the decay rate of the dissociating excited quarkonium state i. These dissociation 
rates for various excited heavy quarkonium states have not been evaluated. It will be inter- 
esting to calculate these rates using explicit quark model wave functions and the formalism 
discussed by Ackleh, Barnes and Swanson pB[ . 

In order to appreciate the effect of dissociation by thermalization on quarkonium survival 
probability, one can consider as an example the case of placing a J/ip in a medium below 
the dissociation temperature of J/ip- If the J/ip is not in thermal equilibrium with the 
medium, then this J/ip system will be stable against spontaneous dissociation. However, 
with the approach of thermal equilibrium (by collision with particles in the medium) at 
a temperature T, the initial J/ip system will evolve into a mixed state with a probability 
distribution to populate different states of the quarkonium system, some of which lie above 
their dissociation thresholds at that temperature. As a result, a fraction of the charmonium 
system will dissociate into open charm mesons, even though the temperature is below the 
J/ip dissociation temperature. Thus, a non-equilibrated J/ip that is stable in the medium 
may become partially unstable against dissociation when it reaches thermal equilibrium. The 
dissociation probability of a quarkonium system depends on its state of thermal equilibrium 
or non-equilibrium. 
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Fig. 6. The fractions of charmonium and bottomium lying above 
the dissociation threshold clS cl function oiT /T^. 
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VII. DISSOCIATION OF HEAVY QUARKONIUM BY COLLISION 



In the hadron thermal bath in which a heavy quarkonium is placed, there will be hadrons 
which will collide with the heavy quarkonium. Dissociation can occur as a result of these 
collisions. 

Previously, we studied the dissociation of heavy quarkonium by collision with light 
hadrons at T = |r5|-[T^. We study in this section the dissociation of J/ip and T by 
collision with pions in a medium at temperature T. The heavy quarkonia under consider- 
ation can be part of a system in thermal equilibrium with the medium. They can also be 
non-equilibrated heavy quarkonia introduced into the medium. The kinetic energy of colli- 
sion would be of the order of the temperature of the hadron gas. As the temperature of the 
medium increases, the quarkonium single-particle energy increases and the energy needed to 
dissociate the heavy quarkonium decreases. If the energy threshold for the dissociation of 
the heavy quarkonium is comparable to the temperature of the hadron gas, the dissociation 
cross section will be large. 

We can calculate the dissociation cross sections of J/ip and T in collision with vr as 
a function of the temperature using the Barnes and Swanson model 0. The calculation 
requires the energies and wave functions of the initial and final meson states, as well as the 
interquark interaction which leads to the dissociation. For the interquark interaction, we 
generalize the temperature-dependent Yukawa and exponential interaction in Eqs. ( ^.81 ) and 
( P.12| ) to the form 

2 ■ 2 \ r 4/i(r) 3m,m/*^ ^I^ttS/^^ J' ^ ' 

For an antiquark, the generator A/2 is replaced by — A-^/2. We have already calculated the 
single-particle energies and wave functions of J / ip and T at various temperatures in Sections 
lT| and 1^ We can calculate the wave functions of vr, open charm, and open bottom mesons 



using the above interaction. The dissociation reactions with the lowest final meson masses 
are J/^^ + vr ^ DD\ D*D, D*D* and T + tt ^ BB*, B*B, B*B*. As the final states involves 
the D* and the B* while the above potential of P7^)- prT^ ) refers to DD or BB at r — »• 00, 
we also need the mass differences of M{D*) — M{D) and M{B*) — M{B), which we shall 
take from the experimental masses at T = 0. According to the spectral analysis by Hatsuda, 
the mass of pion does not change substantially as a function of temperature (see Fig. 3 of 
Ref. [^). We shall keep the mass of pion to be the same as in T = in this calculation. 

The reaction matrix element of Vij is a product of the color matrix element, the flavor 
matrix element, the spin matrix elements, and the spatial matrix elements, as discussed in 
detail in and |T^. With our Gaussian basis states, the spatial matrix elements of Vij 
between the initial and the flnal states can be reduced to the evaluation of Gaussian integrals 
of the Fourier transform of the interaction potential, as shown in Eqs. (48)- (50) of Ref. |]18 . 
For the above Yukawa and exponential interactions, the corresponding integrals are 

/ dq e-(''-''o)V2/3VY.kawa(g) = ^{^2^(3)'!.^, (7.2) 

f dq e-('^-'^")^/2^Ve.ponential(q) = -f^(v^/?)%, (7.3) 
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where 

J„= r r'+^e-'^'''/^-^''' sin qor dr. (7.4) 
Jo 

The integrals X_i and Tq are 

T_i = {exp{zi2}erfc(2;i) - exp{z^}erfc(2;2)} (7.5) 



and 



2^0 = yf^\/f {^iexp{zi}erfc(zi) - ^2 exp{z2}erfc(z2)} , (7.6) 



where 



and 



= (7.8) 
Using these results, the reaction matrix elements can be evaluated and the dissociation 



cross section can be calculated, as in |T^. The sum of dissociation cross sections for vr + 
J /ip ^ DD*, D*D, D*D* are shown in Fig. 7 for different temperatures T/T^, as a function 
of the kinetic energy. Eke = \/s — Ma — Mb where Ma and Mb are the masses of the 
colliding mesons. In Fig. 8 we show similar total dissociation cross sections for vr + T — >■ 
BB* , B*B, B*B* for various temperatures, as a function of Eke- Each cross section curve 
is the average of the results from the "prior" and "post" formalisms, as explained in Refs. 
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Fig. 7. Total dissociation cross section of J /if) in collision with vr 
for various temperatures as a function of the kinetic energy Eke- 

We observe in Fig. 7 that the dissociation cross sections increase as the temperature 
increases. Such an increase arises from decreasing threshold energies as the temperature 
increases. At T/T^ = 0.95 the reaction ir + J/ip D + D* is exothermic and the total 
dissociation cross section diverges at E^e = 0. There is another cross section maximum of 
about 10 mb at Eke ~ 0.17 GeV, which arises from the vr + J ^ D* + D* reaction. As 
the temperature decreases, the cross section decreases. At the temperature of T/Tc = 0.7, 
the maximum cross section decreases to about 5 mb at Eke ~ 0.3 — 0.5 GeV. At T = 0, the 
maximum cross section decreases to about 1 mb at Eke ~ 0.7 GeV and the threshold is at 
0.64 GeV 
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We present similar results for the dissociation of T in collision with vr in Fig. 7b. The 
reactions remain endothermic even for T close to Tc. For T/Tc = 0.95 the maximum cross 
section of about 15 mb occurs at Eke ^ 0.2 GeV. As the temperature decreases, the 
maximum of the cross section decreases and moves to higher kinetic energy. For T/Tc = 0.75, 
the maximum cross section is about 4 mb and is located at Eke = 0.55 GeV. At T = 0, the 
maximum cross section is about 0.6 mb at Eke ~ 1-04 GeV and the threshold is at about 
1.00 GeV []TB[. The dissociation cross section is a sensitive function of the threshold. 

The results in Figs. 7 indicate that over a large range of temperatures below the phase 
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transition temperature, dissociation cross sections of J /if) and T in collisions with vr are 
large. j q 
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c 

Fig. 8. The average dissociation cross section of J/ip and T in 
collision with vr as a function of the temperature. 

In a hadron gas, pions collide with the heavy quarkonium at different energies. We can 
get an idea of the energy-averaged magnitude of the dissociation cross section by treating 
the pions as a Bose-Einstein gas at temperature T. In Figure 8, we show the quantity {av) 
which is the product of the dissociation cross section and the relative velocity, averaged over 
the energies of the pions. The quantity {av) is about 3 mb at T/T^ = 0.70 and rises to 
about 4.5 mb at T/Tc = 0.95 where the value of Tc has been taken to be 0.175 GeV 0. 




VIII. DISSOCIATION OF HEAVY QUARKONIUM IN COLLISION WITH PIONS 

We can estimate the survival probability of a heavy quarkonium in a hot pion gas in the 
presence of this type of collisional dissociation. If we represent the survival probability S by 
exp{— /}, the exponential factor / is given by 

/" ''"freeze 

/=/ {av){T)p{r)dr (8.1) 

J To 

where a is the dissociation cross section, v is the relative velocity between tt and the heavy 
quarkonium, p(r) is the density of vr at the proper time r, and Tq and Tfreeze are the initial 
proper time and the freeze-out proper time respectively. The quantity (av) in Fig. 9 can be 
represented approximately by 

{av) ^ {av),iT/n), (8.2) 
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where {av)c ~ 4.5 mb for J /if) or T. Assuming a Bjorken type of expansion in which the 
density p(r) is proportional to and T{t) is proportional to r~^/^, we obtain 

rr . ( / \ 1/3' 

\ o J 1 / Pfrccz 



/=McPoro(^^j3|l-(^^) }, (8.3) 

where po is the pion density after all the pions are initially formed and the heavy quarkonium 
is thermalized in the medium, Tq is the corresponding temperature, and phceze is the freeze- 
out pion density. Continuing the assumption of Bjorken expansion, the initial pion density 
is given by 

Po = — r-T, (8.4) 
TodyA 

where dN/ dy is the rapidity density of pions and A is the overlapping area in the heavy- ion 
collision, which depend on centrality. We therefore obtain 

To make some estimates on the survival probability of J/ ip and T in a hot hadron gas, we 
can take the quantity tq to be (1 fm/c + 2 fm/c + 2R/'~f), where 1 fm/c is the formation time 
for pions and heavy quarkonium, the additional 2 fm/c is the time for the thermalization of 
hadrons, and ^R/j is the approximate time spread of the first and the last nucleon-nucleon 
collisions in a nucleus-nucleus collision {R is the radius of the colliding nucleus and 7 is the 
relativistic factor in the CM. frame). We can take pfreeze, the freeze-out pion density, to be 
0.5/fm^, corresponding to an average pion separation of about 1.3 fm, and Tq/Tc ~ 1. 

For the most central Pb-Pb collision at 158A GeV, dN^h/dy ~ 450 [g^l, and dN/dy ~ 675. 



Then for this central collision, / = (crt')c(2.40/fm^) and the heavy quarkonium survival 
probability is ^ = e'^-^^ = 0.34. 

For the most central Au-Au colhsion RHIC at ^/sNN = 200 GeV, dNch/dy ~ 650 [|28| , 
and dN/dy ~ 975. Then for this central collision, / = (cr?;)c(7.13/fm^) and the heavy 
quarkonium survival probability is 5* = e~^'^^ = 0.04. There is a substantial absorption of 
both J/tp and T by the hot pion gas at this temperature. 



IX. DISCUSSIONS AND CONCLUSIONS 

We study the dissociation of a heavy quarkonium in high temperatures. The temperature 
of the medium alters the gluon and quark fields and changes the interaction between a 
heavy quark Q and antiquark Q placed in the medium. We first examine the effects of the 
temperature of the medium on the quarkonium single-particle states by using a potential 
inferred from lattice gauge calculations and evaluate the dissociation temperature at which 
the quarkonium dissociate spontaneously. We include spin-dependent interactions and take 
into account the selection rules. We confirm the general features of the results of Digal 
et al. but there are some differences in the dissociation temperatures, which depend on 
the potential and interactions, as well as on the selection rules. We find that the selection 
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rules change the dissociation temperatures substantially for charmonia but only slightly for 
bottomia. 

The results of Digal et al. and our work indicate that most, if not all, of the dissociation 
temperatures of the charmonium states are below T^, but the dissociation temperatures of 
a number of low-lying bottomium states lie above Tc because of the strong binding of these 
states. 

The quarkonium placed in a medium can collide with particles in the medium to reach 
thermal equilibrium with the medium. A heavy-quarkonium in thermal equilibrium can 
dissociate by thermalization as there is a finite probability for the system to be in an excited 
state lying above its dissociation threshold. We find that the fraction of the quarkonium lying 
above the threshold increases with increasing temperatures and is quite large for charmonium 
at T slightly below T^. 

Dissociation of a heavy quarkonium can occur in collision with hadrons at tempera- 
tures below the dissociation temperatures and the phase transition temperature. As the 
temperature increases, the quarkonium single-particle energies increase, and the threshold 
energies for coUisional dissociation decrease. As a consequence, the dissociation cross sec- 
tions increase as the temperature increases, reaching an average of {av) ~ 3.0 mb for the 
dissociation of J /if) and T in collision with pions just below the phase transition temper- 
ature. We have estimated the survival probability of J /if) and T in collision with pions 
in central Pb-Pb collisions at SPS and RHIC energies and found the survival probability 
to be small. Dissociation of J/ip by collisions with pions at high temperatures may be an 
important source of the anomalous J/ip suppression in Pb-Pb collisions observed at CERN 

It has been suggested that collision of J / ip with "comovers" was the source of the anoma- 
lous J/ip suppression in Pb-Pb collisions p]|-P^. There were some uncertainties as the dis- 
sociation cross section of the "comover" with J /ip was not known and the dissociation cross 
sections for different hadrons are very different. The large threshold energy of 0.64 GeV 
for J/ip dissociation in collision with pions raised questions whether the n-J/ip dissociation 
cross section could be large enough to lead to the anomalous J/%1] suppression. The new 
information from lattice gauge calculations of Karsch et al. provides valuable input to 
infer the in-medium dissociation energy of J/ip. The large cross sections for the in-medium 
dissociation in collision with pions obtained here suggest that further microscopic investiga- 
tions of the J /if) dissociation in collision with hadrons at high hadron temperatures will be 
of great interest. 
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